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o 

^ Abstract. 



The Exeter FCRAO CO Galactic Plane Survey consists of ^^CO and ^^CO (J=l-0) 
observations over the galactic plane covering 55° <£<102Wb\ > 1° and 141° <^<195°, - 
3.5° <b<5.5° with a spatial resolution of ~ 45'' and a spectral resolution of ~ 0. 15kms"^ . 
We will present the methodology of a threshold-based cloud and clump determination 
method which retains hierarchical information, then discuss associating sources with 
clouds in the catalogue. Once complete, this catalogue of clouds and clumps will 
encompass the majority of the Northern Galactic Plane, providing knowledge of the 
molecular structure of the galaxy and the starting point for studies of the variation in 
star formation efficiency. In addition, it will allow us to identify clouds that have no or 
fH little star formation taking place inside them, which are often overlooked in the study 

Q^ of the conditions required for star formation to take place. 

o 

, ^, 1. Introduction 

Cn| One of the principle objectives of any millimetre CO observations is to obtain the prop- 

K>^ erties, such as size, velocity dispersion and mass, of the detected molecular structures 

f^ (e.g. [Brunt et al.|2003] ). In data covering large regions in (f,b,v) space the initial identi- 



y^ fication of structures is non-trivial as there will be a range of scales on which emission 

(T^ is found and several independent regions may be detected along a given line-of-sight 

1^ (LOS) due to the structure of the galaxy. Indeed substructure is also likely to be found 

^ on scales down to the spatial resolution of the observations which will probably not 

O break into discrete populations. Another key objective is to be able to associate sources 

of emission at other wavelengths with CO emission structures so as to gain combined 
information about gas and continuum emission, and to pass distance and/or velocity 
information to 2-dimensional data (e.g. |Urquhart"e"t al.„2008a^. Fo r ease of discussion 



we will use nomenclature similar to that of [Williams et al.|pOOOl ) in that we define a 



C^ cloud to be a region of molecular emission which is isolated at the lowest detectible 

level in the data, a clump to be an intermediate object in the heirarchical substructure 
within a cloud and a core to be a top level structure of a similar size to the beam. 

The needs of the two objectives discussed above have often required different de- 
cisions regarding the identification of molecular structures. Obtaining cloud properties 
and/or understanding galactic scale structure (e.g. spiral arms) requires determination 
of low-level strucutures (i.e. clouds) but is less concerned with substructure. Source 
association requires good detection of dense structures (i.e. cores), but is often less 
concerned with intermediate and larger scale structures. Many previous studies have 
biased their structure determination towards the size scale they are interested in. 

1 



Mottram and Brunt 



Given the continuous hierarchical nature of the structure within molecular clouds, 
any method for automatically identifying objects should take this into account in order 
that the results are neither biased towards or away from certain size scales or incomplete 
at a given scale. [Rathbome et"aL] ( |2009| ) used a modified version of the CLUMPFIND 
algorithm ( [Williams et al. 11994 ) with two diff'erent sets of input parameters in order to 
identify both clouds and clumps. However this did not truly track the hierarchy on all 
scales and was limited by the discrete way in which CLUMPFIND steps in temperature 



while searching for structures. [Rosolowsky et al.| ( |2008p developed a cloud decompo- 
sition method which follows the heirarchical structure of a region using dendrograms 
(see e.g. [Houlahan & Scalo|1992[ ) starting from each local maxima down to a common 
threshold until it merges with a neighbour, then considering the merged structures un- 
til a base threshold T^/„ is reached. At each level objects which do not consist of a 
minimum number of pixels Npi^ within a given temperature diff'erence AT of the local 
peak temperature are rejected and, if connected at some lower threshold above Tmin to 
another local maxima are reassigned to that local maxima, similar to the criteria used 
by |BruntetaLl ( |2003l ). 



2. The Exeter-FCRAO CO Galactic Plane Survey 

The Exeter-FCRAO CO Galactic Plane Survey (GPS) consists of ^^CO and ^^CO (J=l-0) 
fully sampled mapping observations obtained with SEQUOIA on the FCRAO with 45'^ 
and 46^' spatial resolution respectively and ~0.15kms"^ spectral resolution over the re- 
gions 55.6°< i <102.5°, \b\^l.25° and 141.5°< i <193°, -3.5°<b<5.5\ This is broken 
up into roughly 4 regions (see Figure[T]) with slightly diff'erent velocity ranges in order 
to follow galactic emission, and will be made public in the near future (Brunt, Heyer, 
Mottram & Douglas, in prep.). The reg ion between 102.5 ° < £ <141.5° is covered by 
the ^^CO Outer Galaxy Survey (OGS, Heyer et al. 1998 |) which was observed with 
QUARRY on the FCRAO. 

These data, though not covering as large an area, provide a factor of 10 improve- 
ment in spatial resolution and a factor of 4 in spectral resolution compared to those of 
Dame et jL] ( |2001| ). In addition, along with the Galactic Ring Survey (GRS, packson 
et al.„20()6p they provide ^^CO (J=l-0) observations for the vast majority of galactic 



plane in the first and second quadrants. 



3. A scheme for heirarchical object detection 

Our method for object detection has been developed from scratch, but has philosophical 
similarities to the approaches of |Brunt et al.| ( |2003] ) and |Rosolowsky et aL] ( |2008| ). While 
the core methodology is general to any 3 dimensional intensity data, it is tailored to data 
in (f,b,v) space with intensities in temperature. 

We begin by thresholding the data at some minimum temperature Tmin, usually 
defined by some integer value m of the median cTrms noise of the data. Each contiguous 
region is checked to see that it passes a series of criteria, which are designed to ensure 
that only self-contained detections are accepted, rather than noise. These are that the 
region contains at least a minimum number of pixels N^/^ (usually calculated from the 
resolution of the data f ^^^^/^^^v). that the local maximum temperature T^nax within the 
region is detected above some diff'erence from the threshold Tmax-Tmin > ST (usually 
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Figure 1. Integrated maps of ^^CO (J=l-0) emission as observed by the 
FCRAO between 55.6°< £ <193° by the Exeter-FCRAO CO GPS and OGS. 
The boxes indicate the different regions covered: The Vulpecular Rift (green, 
55.6°< ^ <68.3°, |/?|^1.25°, -65kms-i<v<65kms-^), Cygnus (blue, 68.2°< 
£ <73.8°, \b\^125\ -85kms-i<v<45kms-i and 73.7°< £ <102.8°, |/?|^1.25°, 
-105kms-i<v<25kms-i), the OGS (red,102.5°< £ <141.5°, -3.5°</?<5.4°, - 
150kms-^<v<40kms-^), the Extended Outer Galaxy (purple, 141.4°< ^ <175°, - 
3.5°<Z7<5.5°, -105kms-^<v<25kms-^) and the Anticentre (cyan, 175°< £ <193°, 
-3.5°<Z?<5.5°, -65kms-i<v<65kms-i). 



an integer number n times the median cr^^^ noise of the data) and that the intensity 
weighted centroid lies within the identified region. These regions are then termed the 
parent objects using tree terminology. 

Following this, each region is examined to determine the minimum threshold be- 
tween Tfnin and the local Tfnax where it divides into two or more separate objects, each 
of which satisfy the criteria discussed above. If this is the case, new objects are created 
which are daughters of the parent region under consideration, which is then removed 
from the list of objects to search and the search moves on to the next region. Alter- 
natively, if the region does not divide to within ST of the local Tmax, the region under 
consideration is the highest level detectable object relating to that local maximum. This 
process is repeated, including checking the daughter objects, until all regions are fully 
searched. 
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While this process is undertaken, a table is constantly updated which contains the 
cloud number and parent(s) and daughter(s) of every detected object. Next, for each 
parent cloud all top level daughters and the pixels assigned to lower objects in the local 
hierarchy within the parent are identified. Running from high to low temperature, each 
pixel is examined to see if any of the surrounding pixels are assigned to one of the top 
level daughters. If not, the next pixel is examined, if true then the pixel is reassigned 
to the top level daughter. This process is repeated until all pixels are assigned to a 
top level daughter. Where more than one top level daughter has a 'claim' on a pixel 
under consideration the daughter with the smallest temperature diff'erence is chosen. 
This process, similar to the CLUMPFIND algorithm, is performed because the density 
structure of structures within molecular clouds is continuous, with dense cores 'owning' 
part of a lower density envelope which is shared with neighbouring cores. This bijection 
paradigm more accurately represents cloud and core properties than simple clipping 
them at a given threshold (see |Rosolowsky et al.|2008[ for a more detailed discussion). 

This top level daughter deconvolution mask is then used to determine the prop- 
erties of all cores, clumps and clouds. For clumps and cores at lower levels within 
the hierarchy, the pixels belonging to them are simply the sum of all pixels from the 
multiple top level daughters they contain. 



4. Applications 

The following sections detail a few of the potential uses and applications of a large 
catalogue of molecular clouds, clumps and cores, aside from exploring the properties 
of the structures in their own right. 

4.1. Source Association 

The association of sources detected in 2D imaging observations with molecular material 
enables velocity (and thus distance) information to be passed from the cloud to the 
sources associated with it. In addition, information about the star formation content and 
emission at other wavelengths of each molecular structure can be obtained. However, 
many of the current methods used to perform this association involve some degree of 
human decision making and hand-crafting (e.g. |Chapin et al.|2008| ), due to the fact that 
there can often be more than one molecular structure along a given LOS. 

Given the scale of both the region observed by the Exeter-FCRAO CO GPS and 
current infrared, sub-mm and mm surveys (e.g. [Churchwell et al.||2009t |Moore et al. 
2005[[Molinari et al.|2010j a more automated approach is required. ]Brunt et aL ~( 2003 ) 
calculated the expectation N^ that a given location will be within 6y of emission with 
peak temperature Tp by seeing the fraction of randomly located test particles that fulfil 
this criteria. ,Kerton & Brunt] (|2003|) then used this expectation to find the most likely 



association between the clouds identified by [Brunt et al.| ( |2003] ) and IRAS point source 
catalogue (PSC, |Beichman et al.|1 988) sources. 

Adopting this approach, we compare results obtained using this method with those 
of two more hand- worked studies in a region around ^=59°. [Chapin et al.| ( |2008) ) used 
an early reduction of our CO data to obtain velocities for their Balloon-borne Large 
Aperture Submillimeter Telescope (BLAST, [Pascale et al.|2008| ) sub-mm point sources. 
The Red MSX Source (RMS) Survey, which is a multi-wavelength search for young 
massive stars through out the galactic plane |Hoare et al.| ( [2005| ); |Urquhart et aL] ( |2008b| ); 
Mottram et al. (2010), performed their own pointed ^^CO observations to their sources 
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Figure 2. The ve locity diffe rence between velociti es associated with BLAST 
( [Chapin et~aL]|20Q8| ) and RMS ( |Urquhart et al.pOQSal source with those assigned 
by associating the sources to molecular structures in a heirarchical cloud catalogue 
using the methodology of Kerton & Brunt (2003 j. 



in this region ( [Urquhart et al.|200'8a| ) in order to obtain kinematic information for their 
sources. The difference between our results using the method discussed above and 
those obtained by |Chapin et al.| ( |2008) ) and |Urquhart et al.| ([2008a) are shown in figure|2] 
These are small and occur mainly due to the fact that we assign the velocity obtained 
from the centroid of the cloud to the source, while the original determinations used the 
CO velocity along the line of sight of the source. The latter approach has been know to 
cause large differences in derived velocity between nearby sources which are known to 
reside in the same complex. For example, ( [Chapin et al.|[20 08 ) find a velocities in the 
range 26-3 Ikms"^ for sources associated with a molecular complex in the Saggitarius 
spiral arm for which they set a common distance. The method of [Kerton & Brunt ( 2003 ) 
automatically assigns the same velocity to sources associated with the same molecular 
structure. 

The ability to associate sources with clouds will allow us to obtain the star for- 
mation efficiency of the clouds and explore how this varies across the galactic plane. 
It will also allow us to identify those clouds which, despite having sufficient density, 
have little or no star formation taking place within them. The relative frequency of 
such clouds will provide important clues as to how quickly star formation begins with 
molecular clouds, and provide important targets for studies of the earliest stages of star 
formation. 



4.2. Distance Assignment 

While kinematic distances can be obtained to all objects identified in the molecular 
data, the uncertainties associated with this approach are often quite large. An alterna- 
tive approach which is possible with a heirarchical catalogue of molecular structures 
is to associate sources with distances obtained by more accurate methods (e.g. maser 
parallax, Reid et al.[[2009) to clouds, clumps or cores within the molecular material, 
then allow this distance to be passed on to all connected parts of the hierarchy above 
the local minimum confusion threshold (LMCT). The LMCT is simply the level in the 
hierarchy at which more than one conflicting distance measurement is available. As im- 
proved distance information becomes available, the distances passed to the molecular 
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structures and the other sources associated with then can be updated as well, without 
having to revisit the CO decomposition itself. This will also ensure that complexes 
have the same distance, which may not be the case if kinematic distances are used, as 
discussed above. 

4.3. Comparison with H i 

From the regions identified in the deconvolution masks, other ({,b,v) data sets, such as 
Hi data from the Canadian GPS (CGPS, [Taylor et aLl|2003| ) can be compared to the 



clouds, clumps and cores identified in the CO data. This would enable the kinematic 
distance ambiguity to be resolved for these structures in the inner galaxy using H i self 
absorption techniques. It will also allow cross-comparison with catalogues of Hi ab- 
sorption and emission features ( [Gibson et al.|2005a fF), so that large statistical samples 
of clouds can be selected to help study the interplay between neutral and molecular 
material in our galaxy. 



5. Summary 

We have a method for identifying heirarchical structure within molecular {-b-v observa- 
tions, which will be applied to the Exeter-FCRAO CO Galactic Plane Survey. The cata- 
logue produced by this process will be associated with sources identified at other wave- 
lengths, which will allow distances from various methods to be passed to the clouds, 
clumps and cores within it. What is more the star formation efficiency of the molecular 
structures will be obtained and variations explored. Other applications beyond explor- 
ing the directly measured properties of clouds within the catalogue include comparison 
with Hi emission and absorption structures, comparison of distance ambiguity resolu- 
tion using extinction and H i self absorption and obtaining the Schmidt- Kennicutt law 
within the galactic plane. 

Acknowledgments. This work was supported by the Science and Technologies 
Research Council of the United Kingdom (STFC) Grant ST/F003277/1 to the Univer- 
sity of Exeter, Marie Curie Re-Integration Grant MIRG-46555, and NSF grant AST 
0838222 to the Five College Radio Astronomy Observatory. J.C.M. is supported by 
a Postdoctoral Research Associate grant from STFC. C.B. is supported by an RCUK 
fellowship at the University of Exeter, UK. The Five College Radio Astronomy Obser- 
vatory was supported by NSF grant AST 0838222. This paper made use of information 
from the Red MSX Source survey database at www.ast.leeds.ac.uk/RMS which was 
constructed with support from the Science and Technology Facilities Council of the 
UK. 



References 

Beichman, C. A., Neugebauer, G., Habing, H. J., Clegg, P. E., & Chester, T. J. (eds.) 1988, 
Infrared astronomical satellite (IRAS) catalogs and atlases, vol. 1: Explanatory supple- 
ment 

Brunt, C. M., Kerton, C. R., & Pomerleau, C. 2003, ApJS, 144, 47 

Chapin, E. L., Ade, R A. R., Bock, J. J., Brunt, C, Devlin, M. J., Dicker, S., Griffin, M., 
Gundersen, J. O., Halpern, M., Hargrave, P. C, Hughes, D. H., Klein, J., Marsden, G, 
Martin, P. G, Mauskopf, P., Netterfield, C. B., Olmi, L., Pascale, E., Patanchon, G, 



Molecular cloud determination in the Northern Galactic Plane 



Rex, M., Scott, D., Semisch, C, Truch, M. D. P., Tucker, C, Tucker, G. S., Viero, M. P., 
& Wiebe, D. V. 2008, ApJ, 681, 428 

Churchwell, E., Babler, B. L., Meade, M. R., Whitney, B. A., Benjamin, R., Indebetouw, R., 
Cyganowski, C, Robitaille, T. R, Povich, M., Watson, C, & Bracker, S. 2009, PASP, 
121,213 

Dame, T. M., Hartmann, D., & Thaddeus, R 2001, ApJ, 547, 792 

Gibson, S. J., Taylor, A. R., Higgs, L. A., Brunt, C. M., & Dewdney, R E. 2005a, ApJ, 626, 195 

— 2005b, ApJ, 626, 214 

Heyer, M. H., Brunt, C., Snell, R. L., Howe, J. E., Schloerb, E R, & Carpenter, J. M. 1998, 
ApJS, 115,241 

Hoare, M. G., Lumsden, S. L., Oudmaijer, R. D., Urquhart, J. S., Busfield, A. L., Sheret, 
T. L., Clarke, A. J., Moore, T. J. T., Allsopp, J., Burton, M. G., Purcell, C. R., Jiang, 
Z., & Wang, M. 2005, in Massive Star Birth: A Crossroads of Astrophysics, edited by 
R. Cesaroni, M. EelH, E. Churchwell, & M. Walmsley, vol. 227 of lAU Symposium, 
370 

Houlahan, R, & Scalo, J. 1992, ApJ, 393, 172 

Jackson, J. M., Rathborne, J. M., Shah, R. Y., Simon, R., Bania, T. M., Clemens, D. P., Cham- 
bers, E. T., Johnson, A. M., Dormody, M., Lavoie, R., & Heyer, M. H. 2006, ApJS, 163, 
145 

Kerton, C. R., & Brunt, C. M. 2003, A&A, 399, 1083 

Molinari, S., Swinyard, B., Bally, J., Barlow, M., Bernard, J. P., Martin, P., Moore, T., & The 
Hi-GAL Consortium 2010, ArXiv e-prints. l(S)Q5 .3317j 

Moore, T. J. T., Shipman, R. E, Plume, R., Hoare, M. G., & Jps International Collaboration 
2005, in Protostars and Planets V, 8370 

Mottram, J. C, Hoare, M. G., Lumsden, S. L., Oudmaijer, R. D., Urquhart, J. S., Meade, M. R., 
Moore, T. J. T., & Stead, J. J. 2010, A&A, 510, A89 

Pascale, E., Ade, R A. R., Bock, J. J., Chapin, E. L., Chung, J., Devlin, M. J., Dicker, S., 
Griffin, M., Gundersen, J. O., Halpern, M., Hargrave, P. C, Hughes, D. H., Klein, J., 
MacTavish, C. J., Marsden, G., Martin, P. G., Martin, T. G., Mauskopf, P., Netterfield, 
C. B., Olmi, L., Patanchon, G., Rex, M., Scott, D., Semisch, C, Thomas, N., Truch, 
M. D. R, Tucker, C, Tucker, G. S., Viero, M. R, & Wiebe, D. V. 2008, ApJ, 681, 400 

Rathborne, J. M., Johnson, A. M., Jackson, J. M., Shah, R. Y., & Simon, R. 2009, ApJS, 182, 
131 

Reid, M. J., Menten, K. M., Zheng, X. W., Brunthaler, A., MoscadelH, L., Xu, Y, Zhang, B., 
Sato, M., Honma, M., Hirota, T., Hachisuka, K., Choi, Y K., Moellenbrock, G. A., & 
Bartkiewicz, A. 2009, ApJ, 700, 137 

Rosolowsky, E. W., Pineda, J. E., Kauffmann, J., & Goodman, A. A. 2008, ApJ, 679, 1338 

Taylor, A. R., Gibson, S. J., Peracaula, M., Martin, P. G., Landecker, T. L., Brunt, C. M., 
Dewdney, P. E., Dougherty, S. M., Gray, A. D., Higgs, L. A., Kerton, C. R., Knee, 
L. B. G., Kothes, R., Purton, C. R., Uyaniker, B., Wallace, B. J., WilHs, A. G., & 
Durand, D. 2003, AJ, 125, 3145 

Urquhart, J. S., Busfield, A. L., Hoare, M. G., Lumsden, S. L., Oudmaijer, R. D., Moore, T. J. T., 
Gibb, A. G., Purcell, C. R., Burton, M. G., Marechal, L. J. L., Jiang, Z., & Wang, M. 
2008a, A&A, 487, 253 

Urquhart, J. S., Hoare, M. G., Lumsden, S. L., Oudmaijer, R. D., & Moore, T. J. T. 2008b, in 
Massive Star Eormation: Observations Confront Theory, edited by H. Beuther, H. Linz, 
& T. Henning, vol. 387 of Astronomical Society of the Pacific Conference Series, 381 

WilHams, J. P., Blitz, L., & McKee, C. E. 2000, in Protostars and Planets IV, edited by V. Man- 
nings, A. R Boss, & S. S. Russell, 97 

Williams, J. R, de Geus, E. J., & Blitz, L. 1994, ApJ, 428, 693 



